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1. Introduction 

The total cross sections for the processes pp' — HH' + X, where pp' G {qc},c}C},gg,gc},gq} 
and H,H' are heavy particles in arbitrary representations R,R' of the SU{3) colour algebra (top 
quarks, squarks, gluinos...), contain classes of contributions which are enhanced near the partonic 
threshold region j3 = y^l — 4M^/s — 0, with M = {Mh +M///)/2 the average mass of the two heavy 
particles and s the partonic centre-of-mass energy. These corrections arise from soft-gluon radiation 
off initial- and final-state particles (~ a"ln™j8, "threshold logaritms") and exchange of Coulomb 
gluons between the two non-relativistic heavy particles (~ (ajp)", "Coulomb singularities"), and 
should be resummed to all orders in if the hadronic cross section is dominated by the partonic 
threshold region. Resummation of threshold logarithms in the Mellin-moment space formalism 
has been discussed, for example, in |Jl|, |^, and studies of the effects of Coulomb resummation in 
production of top quarks and SUSY particles have been presented in [^, 0, ^. However, the issue 
of factorisation and simultaneous resummation of soft and Coulomb effects in pair production has 
been addressed rigorously only recently in employing effective-theory techniques similar to 
ones previously applied in the context of DIS and Drell-Yan [^, ||]. Here we present the application 
of results obtained in [^] to the process of squark-antisquark production at the LHC, and discuss 
the numerical relevance of soft and Coulomb resummation for theoretical predictions of the total 
cross section (see also [Q]). 



2. Factorisation and RG evolution equations 

In ^ it has been shown that near the partonic threshold s ~ 4M^ the total partonic cross 
section for the pair-production process factorises, at leading-order in j3 and for particles produced 
in an S'-wave state, according to 

V(^'At/) =I^n'(^,At/) / dco Y^Jj,^{E-^)w;^"{(0,^f), (2.1) 

with E = ^Jl — 2M. Here Hiii{M,jjLf ) denotes a process-dependent hard coefficient encoding 
short-distance effects, decomposed over a basis of colour-state operators {cj'j}}- The soft func- 
tion W^i" {(O ^ pif) is defined in terms of matrix elements of soft Wilson lines, decomposed over 
the basis and projected onto the irreducible colour representations Ra of the HH' system, 

R®R' = Y^a^a- The function Jr^ describes the internal evolution of the heavy-particle pair HH' , 
which is driven by Coulomb-gluon exchange. Near the partonic threshold the basis } can be 
chosen such that the soft function W^"{(0,Hf) is diagonal in colour space to all orders in a,, and 



the colour structure of ( |2.1D simplifies significantly [g]. 

The resummation of threshold logarithms follows the formalism of [^, and is a generalisation 
of the Drell-Yan case [S] to arbitrary colour representations /?„ of the heavy-particle pair system. 



From results on the IR structure of general massive amplitudes in QCD 01O| ] it follows that the 
diagonal elements of the hard function //,(M,/ij) = Hii{M,Hf) satisfy a renormalisation-group 
evolution equation given by 

-Hi{M,pLf) = ( 7cusp(C, + C,-0 In f ^ I + 2/ I , (2.2) 
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Figure 1: Schematic representation of the resummation of threshold logarithms via RGBs. 



where /cusp denotes the cusp anomalous dimension, C, and Cr' are the Casimir invariants of the 
initial-state representations, and yy can be written in terms of single-particle anomalous dimensions 
yf = Y + Y + • Froni the scale-invariance of the total hadronic cross section, one can derive 
an analogous evolution equation for the soft function W^'^"(G),/Xy) [^], 

-^W,^«(«,At/) = -2 [(C, + Q)7cuspln +2T^«.j W,'^«(«,At/) 

dva.\lf L VM// 'J 

-2(c,+c.)r„.,r^.- "''°'"'^^'-";'°°"°-''^' . (2.3, 

JO w — w 

with the anomalous dimension given by = ■J^^ + y*^ '' + 7^ '"', and t^'*" defined by the large-^c 
Umit of the Altarelli-Parisi splitting function 



(2.4) 



3. Momentum-space resummation 

In the approach of the resummation of threshold logarithms is obtained directly in momen- 
tum space by calculating the hard and soft functions at scales pLf, and jig respectively, and evolving 



them to the common factorisation scale /X/ using the RG equations ( [2.2| ) and ( [2.3[ ), as schematically 
depicted in figure |l]. The hard scale jih and soft scale ji^ are chosen to minimise radiative correc- 
tions to Hi{M,iih) and W^" {(O , }is) ■ As anticipated above, the solutions of equations CLl) and (2.3) 
are generalisations of the Drell-Yan case discussed in ^ and read 



2 \ -2ar(AiA,M/) 



Hi{M,IXh); 



1 (CO 



W^^''\co,^f) = exp[-45(At,s-,M/)+24':,(At.v,M/)]^f"(5r„Ms);^ [jj ^(^)r(2r]) 



2rj 



,(3.1) 



where 



ariv,iJ.) = -{Cr + Cr') da 



2j8(a,) 



af (v,At) 



«.(v) i8(a,) 
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and5f"(p,/Xi) is the Laplace-transform of W^-'^"( ft), /Xi) with respect to the variables= 

poo 

5f«(p,Ats)=/ d(Oe-'''wl'"i(0,^,,). (3.3) 

J0_ 



The one-loop expression for sf"(5f,,/is) for arbitrary colour representations was computed in 



Notice that the resummed hard and soft functions in eq. (3.1) are formally independent of the 
choice of /i/j and /i.,, but when //(M,/X/,) and sf°{dri,IJ.s) are truncated at a finite perturbative order 
a dependence on the hard and soft scales is introduced, which is however of higher order in a^. 

The resummation of the velocity-enhanced terms (a^/jS)" associated with Coulomb-gluon 
exchange is obtained by relating the potential function Jr^ to the zero-distance Coulomb Green 



function, Jr^ = 2ImGl?l {0,0;E). Using the representation provided in []11[], this gives 



2n 



Im 



2m, 



red 



■ + a,(-DRj 



2% My' J 



1, , A 0Cs{-DrJ Y 

:z + Ye + ¥\ 1 = 

2 2V-£/(2mred)yJ 



(3.4) 



with Dr^ = {Cr^ -Cr-Cr')/2 and mred = MeMh' / {Mh + Mh') ■ The second term on the first hne 



of (3.4) encodes the contribution to the cross section from a single Coulomb-exchange diagram, 
while the second line accounts for the contribution of two or more gluons. Beside modifying the 
shape of the cross section above threshold {E > 0), the resummation of Coulomb effects leads to 
the appearance of heavy-particle bound states below threshold for an attractive Coulomb potential 
{Dr^ < 0). In the limit of vanishing width, Tu = Tfj' = 0, and for £ < 0, Jr^ reduces to 



JrAE) = 2Y^ 



;i=l 



mga,{-DR^] 
2n 



5\E + ma 



2n 




-Dr.)- 



(3.5) 



4. Threshold resummation for squark-antisquark production 

We now apply the results presented in the previous section to the NLL soft resummation and 
Coulomb resummation of the squark-antisquark total production cross section. The resummed 
cross section for the partonic channel pp' is 



l^Hf^\m^,Hf) I dco l^jR^{E-'l-)W^''^^'^{co,Hf). (4.1) 



The NLL resummed hard and soft functions in ( |4.1[ ) are given by the expressions in eq. ([3.1|), where 
the cusp anomalous dimension /cusp is included at the two-loop accuracy, the anomalous dimen- 
sions yj , at one-loop accuracy, and Hi{M,}ik) and ^"{d-q^fXs) are replaced by their tree-level 
expressions Hf'\M ^jih) and sf"^^\drj,IJ.s) = 1 (for more details on the reorganised perturbative 
expansion in presence of soft and Coulomb resummation we refer to [^]) . The resummed cross 
section ( p~T| ) is then matched onto the full fixed-order NLO result 



(4.2) 
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using the parametrisation of the NLO cross section extracted in [12] from numerical results ob- 
tained with Prospino [13]. Our final prediction for the hadronic cross section is given by the con- 
volution of (gj) with the MSTW2008 NLO PDFs. The default value for the soft scale, denoted 
with fis, is chosen such that the one-loop soft correction to the hadronic cross section is minimised. 
This is analogous to the procedure introduced in and yields jig = 123 — 455GeV for the squark 
mass range = 200 GeV — 2TeV. Furthermore, we evaluate the running coupling in Jr^ at the 
scale He = Max{2m^j8,Cirm^o:i(/Xc)}- This choice is motivated by the characteristic virtuality 
of a Coulomb gluon, ~ "J^i^^, and by the inverse Bohr radius of the first HH' bound state, 
l/rg = Cpniqas/l. Finally, we identify the hard and factorisation scale, /i/, = /ij, and choose the 
latter of order niq. This means that the resummed hard function in ( |4.1[ ) reduces to its tree-level ex- 
pression, since the functions vanish for this particular scale choice. The function Hj^\mq, Hf) 



(0) 




"-[GeV] 



Figure 2: Corrections to the NLO squark-antisquark 
cross section from soft-Coulomb resummation as func- 
tions of the squark mass niq. 



is extracted from the expansion at threshold of the tree-level cross section a, — 

The numerical effects of soft-Coulomb 
resummation on the total hadronic cross 
section, normalised to the full NLO result, 
are shown in figure ^, for y/s = 14 TeV, 
degenerate squark masses, mg/mq = 1.25 
and jtt/ = mq. The four curves represent 
different contributions contained in the re- 
summed cross section (pTII). The low- 
est curve ("EFTsoft", dashed blue) shows 
the effect of NLL soft-gluon resummation 
without any Coulomb effect, correspond- 
ing to setting Jr^ to its tree-level value 
in "EFTsoft+c" (dot-dashed red) in- 

cludes, on top of soft effects. Coulomb 
resummation, but no interference of the 

two. This is partially included in the third curve ("EFTnll"> solid black) convoluting the one-loop 
Coulomb contribution (second term in (3^)) with the full NLL soft function (see also [^). Finally, 
the contribution to the cross section from bound states below threshold is added in the last curve 
("EFTnll+bs"> dotted magenta). The pure soft resummation is in qualitative agreement with the 
Mellin-space results quoted in [^], while the choice of a running scale in the Coulomb corrections 
leads to much larger effects than in [^, where at the scale pLf was used. Also, the soft-Coulomb 
interference and bound-state contributions shown in figure ^ turn out to be comparable in size to 
pure soft and Coulomb effects, and are therefore necessary for precise theoretical predictions. 

The reduced scale dependence of the resummed cross section is presented in figure ||. The 
left plot shows the LO, NLO and resummed result as functions of the factorisation scale /X/, for 
mq = \ TeV and pLf, = /ly . The green band represents the uncertainty associated with the choice of 
the soft scale, which varies in the interval /l,v/2 < pL^ < ^P^s- The scale dependence is clearly reduced 
for large jJ./, but only mildly improved at small values of the scale, where ln(/ii//X/) = 0(1) and 
the effect of soft resummation is negligible. It can be argued that in this region (/X/ < 0.2m^) the 
identification of the hard scale with the factorisation scale is not justified, and the two scales should 
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Figure 3: Scale dependence of the LO (solid red), NLO (dotted blue) and resummed cross section (green 
band) for different choices of jJ-s and jlf,- See the text for explanation. 

be kept separate. This is shown in the right plot in figure |3[ In this case the green band is obtained 
by varying /X/, and {Xg, independently from the factorisation scale, in the intervals < /X/, < 4mq 
and fLs/2 < /i.v < 2/I.v, respectively. While the cross section is not significantly modified with respect 
to the left plot at large scales, in the region /ij ^ the appearance of large logarithms of /X/,//^/, 
and their resummation in Hf^^^{m^,jj.f), leads to sizeable corrections to the cross section. 
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